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We present new results of the search for WH — » ivbb production in pp collisions at a center of mass 
energy of y/s = 1.96TeV, based on a datasct with integrated luminosity of 0.44 fb~^. Wo combine 
these new results with previously published searches by the DO collaboration, for W H and ZH 
production analyzed in the $j,hb final state, for ZH {-^ £^£~bb) production, for WH{^ WWW) 
production, and for H {—^ WW) direct production. No signal-like excess is observed either in the 
WH analysis or in the combination of all DO Higgs boson analyses. We set 95% CL. (expected) 
upper limits on cr(pp — » WH) x B{H — » 66) ranging from 1.6 (2.2) pb to 1.9 (3.3) pb for Higgs 
boson masses between 105 and 145 GeV, to be compared to the theoretical prediction of 0.13 pb for 
a standard model (SM) Higgs boson with mass niH = 115 GeV. After combination with the other 
DO Higgs boson searches, we obtain for mn = 115 GeV an observed (expected) limit 8.5 (12.1) 
times higher than the SM predicted Higgs boson production cross section. For mn = 160 GeV, the 
corresponding observed (expected) ratio is 10.2 (9.0). 

PACS numbers: 13.85Qk,13.85.Rm 



Spontaneous electroweak symmetry breaking in 
the standard model (SM) provides an explanation 
for the masses of the elementary particles, other- 
wise massless in the unbroken gauge theory. Its 
success, in particular in explaining the mass of the 
electroweak vector bosons, awaits one last but nec- 
essary experimental confirmation: the observation 
of the Higgs boson, which is a scalar particle asso- 



ciated with the symmetry breaking. For Higgs bo- 
son searches, the most sensitive production chan- 
nel at the Tcvatron for a Higgs boson with mass 
below 130 GeV is the associated production of a 
Higgs boson with a W boson. All possible chan- 
nels, however, must be studied to gain sensitivity 
through their combination. 

At a center-of-mass energy of y/s = 1.96 TeV, 
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three pp WH searches have already been pub- 
hshed or submitted for pubhcation, one [l[ using a 
subsample (0.17 fb~^) of the dataset used in this 
letter, while the two others are from the CDF col- 
laboration: one uses 0.32 fb~^ 0] of data, the other 
updates it using improved analysis techniques and 
a larger dataset based on 1.0 fb^^ of integrated 
luminosity 

For this WH analysis we require one high trans- 
verse momentum (pr) lepton (e or fi), missing 
transverse energy l^rp to account for the neutrino 
in the W boson decay, and exactly two jets with 
at least one of them being identified as originat- 
ing from a bottom (b) quark jet ( "^-tagged" ) , as 
detailed below. The dominant backgrounds to 
WH production are W+ heavy-flavor production, 
top quark pair production (tt), and single top 
quark production. This analysis uses a dataset of 
0.44 fb~^. Compared to the previous DO result, 
the 6-jet identification has been optimized, and the 
muon channel has been added. 

The result of this search is then combined with 
previously published searches by the DO collabo- 
ration with a similar luminosity. These searches 
cover WH and ZH production analyzed in the 
^Tbb final state [J, ZH i+e-bb) produc- 
tion d, WH{^ WW+W ) production and 
H W+W-) direct production 0. In the fol- 
lowing, the particle charges will not be mentioned 
explicitly, except when needed to resolve potential 
ambiguity. We first describe the WH analysis in 
detail, then the full combination of results. 

The WH analysis relies on the following compo- 
nents of the DO detector 

i) a. central-tracking system, which consists of a sil- 
icon microstrip tracker (SMT) and a central fiber 
tracker, both located within a 2 T superconduct- 
ing solenoidal magnet; 

a) SL liquid- argon/uranium calorimeter with a cen- 
tral section (CC) covering pseudorapidity [l3| 
< 1.1, and two end calorimeters (EC) extend- 
ing coverage to |r/| ~ 3.2, all housed in separate 
cryostats, and with scintillators between the CC 
and EC cryostats providing sampling of develop- 
ing showers at 1.1 < \rj\ < 1.4; 
Hi) a muon system, which surrounds the calorime- 
ter and consists of a layer of tracking detectors 
and scintillation trigger counters before 1.8 T iron 
toroids, followed by two more similar layers behind 



the toroids. 

We reject data periods in which the quality of 
the data in the tracking, the calorimeter, or the 
muon system is compromised. The luminosity is 
measured using plastic scintillator arrays located 
in front of the EC cryostats, covering 2.7 < \r]\ < 
4.4. The uncertainty on the measured luminosity 
is 6.1%. The -|- jets candidate events must pass 
one of the triggers which require, for the e channel, 
at least one electromagnetic (EM) object, and for 
the fj, channel, at least one muon object or a trigger 
requiring a muon and a jet in the final state. 

The event selection for the WH analysis requires 
one lepton candidate with transverse momentum 
Pt > 20 GeV, > 25 GeV, and exactly two jets 
with Pt > 20 GeV and \ri\ < 2.5. Only events 
having a primary z vertex within ± 60 cm of the 
nominal interaction point are accepted. If the lep- 
ton is an electron, it is required to have \ri\ < 1.1. 
If it is a muon the requirement is \r]\ < 2.0. 

Electrons are identified in two steps. The pre- 
selected electron candidates (seeded by an energy 
cluster in the EM calorimeter) are first required 
to satisfy identification (ID) criteria: (a) a large 
fraction of their energy deposited in EM layers, 
i.e. EMF > 0.9, (b) low fractional energy de- 
posited around the expected electron energy de- 
position, and (c) spatial energy distribution in the 
EM calorimeter consistent with that of an elec- 
tron. These criteria define "loose" electrons. The 
loose electrons are then tested with a likelihood 
algorithm, optimized on Z — > ee samples, and 
which takes as input seven quantities sensitive to 
the EM nature of the particles ll|. If they sat- 



isfy the likelihood requirement, they are accepted 
as final ("tight") electrons for the analysis. The 
efficiencies of the ID and likelihood requirements 
are determined from a dielectron sample in which 
we select a pure set of Z events. The combined 
reconstruction and ID efhcicncy is found to be 
(95.4 ± 0.4)%. The hkelihood efficiency for elec- 
trons is (92.0 ±0.3)%. 

Muons are reconstructed using information from 
the muon detector and the central tracker. They 
are required to have hits in all layers of the muon 
system inside and outside the toroid. The superior 
spatial resolution of the central tracker, inside the 
strong solenoidal magnetic field, is used to improve 
the accuracy of kinematic properties of the muon 
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and to confirm that the muon originated from the 
primary vertex. A veto against cosmic-ray muons 
based on the timing of hits in the muon-system 
scintillator detectors is applied. Quality criteria 
on the associated central track are also applied 
to reject the majority of background muons: a 
small track impact parameter (dca) compared to 
its resolution (adca) is required, dca < 3adca, to 
reject muons originating from semi-leptonic decays 
of heavy-flavor hadrons which constitute the main 
background. Such background muons have a lower 
transverse momentum spectrum and are not typi- 
cally isolated due to jet fragmentation. A loose iso- 
lation criterion is defined using the spatial separa- 
tion AR — y%A^p + (A^)2 between a muon and 
the closest jet in the r]-(p plane, where ip is the az- 
imuthal angle, we require AR > 0.5. Tighter muon 
isolation criteria are defined by requiring that the 
scalar sum of the transverse energy of calorime- 
ter clusters in a hollow cone (0.1 < AR < 0.4) 
around the muon divided by the px of the muon 
be less than 0.08, and the scalar sum of the trans- 
verse momenta of all tracks within a cone of radius 
AR = 0.5 around the muon divided by the px of 
the muon be less than 0.06. The track matched to 
the muon is excluded from this sum. 

The jets are reconstructed using a cone algo- 
rithm 18] with a radius of AR = 0.5. We apply 
standard DO jet-ID criteria to avoid fake jets which 
occasionally originate from noise in the calorime- 
ter, i.e., the energy fraction in the EM layers of a 
jet is required to be 0.05 < EMF < 0.95 and the 
energy fraction in the CH section of the calorimeter 
is required to be < 0.4. The difference in efficiency 
of the jet-ID requirements between data and simu- 
lation is quantified in the overall jet reconstruction 
efficiency scale factor to which a systematic uncer- 
tainty of 5% (per jet) is assigned. 

The multijet background is estimated from the 
loose and tight e or /i final samples, as described in 
Ref. pjj] using the following probabilities. We de- 
termine from the data the probability Pi^osci^'ight 
for a "loose" lepton originating from a jet to pass 
the tight lepton requirements. This is done sepa- 
rately for the electron and the muon channel and 
this probability is determined as a function of the 
Pt of the candidate lepton. The sample of mul- 
tijet events containing a loose lepton is selected 
with kinematic criteria that ensure negligible con- 



tamination of real leptons. We also determine the 
same type of probability PioosoLtight ^'^^ * genuine 
isolated lepton from Z £~^£~ samples. With 
these two probabilities and the numbers of loose 
and tight W + 2 ^et candidates, we determine the 
number of multijet background events in our sam- 
ple, bin-by-bin, for every differential distribution. 

To select W boson decays, we require 1^t> 
25 GeV. The is calculated from the calorime- 
ter cells except for unclustered cells in the out- 
ermost layer of the calorimeter (coarse hadronic 
layer, CH) and is corrected when one or several 
muons are present. All energy corrections to elec- 
trons or jets are also propagated into the l^^- The 
transverse mass of the W boson candidates in the 
W+ jets sample is reconstructed from the lepton 
and missing transverse energies. Its distribution is 
shown in Fig. [1] and compared with the sum of con- 
tributions from multijet events with misidentified 
leptons and from SM processes which are obtained 
from simulated events. 
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FIG. 1: Distributions of the transverse W boson mass 
compared to the simulated expectation in the W + 2 
jet event sample. The simulation is normalized to the 
integrated luminosity of the data sample using the ex- 
pected cross sections taking into account all the other 
backgrounds (the fraction of WH events is negligible 
before 6-tagging). 
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The following processes are simulated with the 
PYTHIA MC event generator version 6.202, 
making use of the CTEQ5L [l3[ leading-order par- 
ton distribution functions: inclusive production of 
W ^ e/^J./T + Z ^ ee/^i^/rr; WW, WZ, ZZ; 
tt e/fi/r + jets production (lepton-|-jets and 
dilepton channels), WH e/ ^/t + v + hh produc- 
tion. The single top quark processes are generated 
using COMPHEP [lj|. 

Throughout this Letter, 'W-|-jets" simulated 
events refer to events with a W produced in as- 
sociation with light-flavor jets (originating from u, 
d, s quarks or gluons; generically denoted by j) or 
charm jets (originating from a c quark). They con- 
stitute the dominant background before 5-tagging. 
and are generated with ALPGEN [I^ (interfaced to 
PYTHIA for showering and fragmentation), since 
ALPGEN has a more complete simulation of pro- 
cesses with high jet multiplicities. The generation 
is based on W+ 2 jets {Wjj) processes, includ- 
ing the charm quark (c) processes Wcc and Wcj. 
The Wbb events are generated separately requiring 
two b parton jets with pt > S GeV separated by 
Ai? > 0.4; its NLO cross section is obtained using 
MCFM [lit . 

These simulated backgrounds are absolutely 
normalized (according to NLO cross sections) with 
the exception of the W+ jets sample which is nor- 
malized to the data after subtraction of all the 
other backgrounds. The systematic uncertainty 
on the NLO cross sections of these processes is 6- 
18%, depending on the process. All these events 
are processed through the DO detector simulation, 
based on GEANT [17|, and the reconstruction soft- 
ware. The simulated events are then weighted by 
the trigger efflciency and by the data/simulation 
ratio of all the selection efficiencies. The shape 
of the distribution of the transverse mass of the 
W candidates (Fig. [T]) is well reproduced by the 
simulation of the W + jets processes, after adding 
the multijet background and the other SM back- 
grounds. 

To identify heavy-flavor jets we use a ^-tagging 
algorithm which computes a probability correlated 
to the b quark lifetime [l9l |. The requirements on 
the "jet lifetime probability" (JLIP) have been op- 
timized for events with one or two 5-jet candidates 
by maximizing the sensitivity to the Higgs boson 
signal. The requirement is first set to 1%; if two 



jets are tagged the event is selected as double b- 
tagged (DT). Otherwise the requirement is tight- 
ened to 0.1% and if exactly one jet is tagged the 
event is selected as single 6-tagged (ST). In this 
way the single and double fe-tagged subsamples are 
independent, which simplifies their combination. 
The mistag rate (tagging of light flavor jets) ob- 
tained in these samples are approximately equal 
to the corresponding JLIP requirements, while the 
efficiency for correctly identifying a genuine & jet 
("^-tagging efficiency") is (55±4)% and (33±4)%, 
respectively. These efficiencies were determined 
with central "taggable" jets (I77I < 1.2) having a 
transverse momentum of 35 < < 55 GeV. A jet 
is "taggable" if at least 2 tracks (one with pT > l 
GeV, the other with pr > 0.5 GeV) and > 1 SMT 
hits are inside the Ai? < 0.5 cone defining the jet. 
The jet taggability is typically 80% in a two-jet 
sample with an uncertainty of 3%. 

For each tagged jet in the simulation, we apply 
the ratio between the expected taggability times 
6-tagging efficiency in data and in simulation to 
reweight the simulated events. For the tagging ef- 
ficiency of simulated 5 or c jets, we use pr — 
dependent data vs. simulation scale factors, de- 
termined from real b jets |19i] . In the simulation, 
the tagged light flavor jets are weighted to repro- 
duce the mistag rate as measured in data using 
dedicated samples fioj . 

With the above selection criteria, we observe 
137 W + 2 jet events having exactly one 6-tagged 
jet (ST sample) and 30 events having both jets b- 
tagged (DT sample) . In these samples the multijet 
background is estimated using as a loose sample 
the W + 2iet ST (DT) sample in which the lepton 
is selected using the loose lepton-ID criteria. The 
distribution of the invariant dijet mass of 14^ -I- 2 
jet events for the ST and DT samples is shown in 
Fig. and b. The data are compared to the sum 
of the simulated SM processes added to the multi- 
jet background. The agreement indicates that the 
simulation describes the data well. 

The different components of the background are 
shown in Table HI The small expected contribu- 
tions from a 115 GeV Higgs are also shown, but 
no excess above the standard model backgrounds 
is visible in these distributions, so we proceed to 
set limits from these distributions, after systematic 
uncertainty evaluation. 
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The experimental systematic uncertainties on 
the efficiencies and those due to to the propaga- 
tion of other systematic uncertainties (trigger, en- 
ergy caUbration, detector response) which affect 
the signal and SM backgrounds are the following 
(ranges indicate different values for the e and ^ 
channel): (2-3)% uncertainty from the trigger effi- 
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FIG. 2: Dijet mass distributions for the W + 2 jet 
events (a) when exactly one jet is tightly 6-tagged 
and (b) when the two jets are loosely 6-tagged (see 
text). The data are compared to Wbb, ti, VK-|-jets and 
other smaller expectations. The background labeled as 
"other" in the figure is dominated by single top quark 
production. 
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WH 
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0.49 


± 


0.07 


0.43 
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U.UD 


WW,WZ,ZZ 148.7 ± 23.8 


5.3 


± 


0.8 


2.0 


± 


0.4 


Wbb 


116.3 ± 18.6 


22.3 


± 


4.4 


14.4 


± 


3.2 


tt 


87.6 ± 8.6 


21.0 


± 


4.3 


12.6 


± 


2.7 


Single top 


41.2 ± 5.3 


10.0 


± 


4.8 


3.7 


± 


0.7 


Multijet 


984 ± 153 


22.8 


± 


7.5 


1.5 


± 


0.6 


W/Z + jets 


6908 ± 1076 


57.7 


± 


10.3 


4.1 


± 


0.7 


Total expect. 


8286 


139.6 


± 


28.5 


38.7 


± 


5.8 


Observed Ev. 


8286 


137 






30 







TABLE I; Summary table for the £ (e and fi) + 2 jets -|- 
final state. Observed events in data are compared 
to the expected number ofW + 2 jet events before and 
after 6-tagging in the simulated samples of WH, di- 
bosons, Wbb production, top production {ti and single 
top), multijet background, and "W /Z+ jets" produc- 
tion. In the pre-tagged sample the W /Z + jets contri- 
bution is normalized such that the total expectation is 
normalized to the data. 



ciency, (3-4)% uncertainty for the lepton identifi- 
cation and reconstruction efficiency, (3-4)% for the 
lepton energy scale and resolution, 5% for the jet 
identification and reconstruction efficiency, 5% for 
the modeling uncertainty of the jet multiplicity in 
the simulation, (5-12)% due to the jet energy cali- 
bration uncertainty, 3% for the jet taggability, and 
(5-6)% for the 5-tagging efficiency; for the light 
quark jets these uncertainties are 9% (DT) and 
13% (ST). In summary, for WH production and 
simulated backgrounds, the experimental system- 
atic uncertainty is (16-19)%. The multijet back- 
ground, determined from data, has an uncertainty 
of 25%. The systematic uncertainty on the cross 
section of the simulated backgrounds is 6-18%, de- 
pending on the process. The uncertainty on the 
luminosity is 6.1%. 

The limits for WH production are obtained us- 
ing the CLs method [20|, [21 1 taking the dijet in- 
variant mass of the bb system as the final dis- 
criminating variable. It is performed on the ST 
and DT samples of the e and ^ channels inde- 
pendently (four analyses), which are then com- 
bined. The CLs approach is based on the like- 
lihood ratio test statistic, Q = L{s + b)/L{b) ~ 

e~'° + '''(s+b)" /e-^b)" 



-/- 



where s and b arc the ex- 
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FIG. 3: LLR distributions obtained with the CLs 
method for the combination of the ST and DT samples 
in the WH channel. 



pected numbers of signal and background events 
while n is the number of data events. For com- 
putational ease, the log-likelihood ratio LLR{n) = 
—2 ln{Q) is used. In order to exploit the shape in- 
formation of the final discriminating variable, as 
well as combine the different channels, the LLR 
values per bin and for all channels are added. Sys- 
tematic uncertainties are incorporated into the sig- 
nal and background expectation using Gaussian 
sampling of individual uncertainties. Correlations 
between uncertainties across channels are handled 
by varying simultaneously the fluctuations of iden- 
tical sources of all channels. The 95% C.L. lim- 
its are determined by raising the signal cross sec- 
tions until the ratio of probabilities for the sig- 
nal+background hypothesis to the background- 
only hypothesis falls below 5%. 

Figure [3] shows the LLR distributions for 
the WH combined result. The LLR values 
for the signal -|-background hypothesis {LLRg+b), 
background-only hypothesis {LLRb), and the ob- 
served data [LLRobs) are shown. The quantities 
LLRs+bi LLRb, and LLRobs are obtained by set- 
ting n = s -\- b OT n (observed) into LLR{n). 
The shaded bands represent the one and two stan- 
dard deviation (a) departures for LLRb- These 
distributions can be interpreted as follows: The 
separation between LLRb and LLRs+b provides a 
measure of the discriminating power of the search; 
the width of the LLRb distribution provides an es- 
timate of the sensitivity of the analysis to a signal- 
plus-background-like fluctuation in data, taking 
account of the systematic uncertainties; the value 



of LLRobs relative to LLRs+b and LLRb indicates 
whether the data distribution appears to be more 
signal-like or background-like, and the significance 
of any departures of LLRobs from LLRb can be 
evaluated by the width of the LLRb distribution. 

The observed (expected) combined upper limits 
obtained at 95% C.L. on a(pp WH) x B{H 
bb) range from 1.6 pb to 1.9 pb (2.2 pb to 3.3 pb) 
for Higgs boson masses between 105 and 145 GeV 
and are displayed in Fig. |4l They are also given 
in Table El together with the ST and DT subchan- 
nel limits and the ratios of all these limits to the 
predicted SM cross section. These new WH up- 
per limits are compared in Fig. S] to the previ- 
ously published results on WH production from 
DO on 0.17 fb~^ of data in the electron channel 
only [H and CDF (0.32 fb"! e+ii channels) 0. 
The improvement in sensitivity obtained with this 
analysis is clearly visible in the region where the 
Tevatron is most sensitive to a Higgs boson with 
mass in the 115-135 GeV range. The result is also 
compared to the CDF result recently submitted 
for publication on 1.0 fb~^ of data [3|, showing 
comparable expected sensitivity when taking into 



WH 1 Higgs mass [GeV] 


105 


115 


125 


135 


145 


ST observed a x B 
ST expected a x B 


6.62 
8.11 


5.74 5.17 
6.94 6.10 


4.79 
4.90 


4.74 
5.08 


DT observed a x B 
DT expected a x B 


2.21 
3.55 


2.12 
3.07 


2.25 
2.89 


1.98 
2.43 


1.97 
2.58 


ST4-DT observed a x B 
ST-(-DT expected a x B 


1.92 
3.25 


1.71 
2.83 


1.79 
2.53 


1.64 
2.16 


1.77 
2.21 


ST observed ratio to SM 
ST expected ratio to SM 


34.9 
42.8 


44.9 
54.4 


65.5 
77.3 


112.8 
115.4 


250.7 
268.5 


DT observed ratio to SM 
DT expected ratio to SM 


11.7 
18.8 


16.6 
24.1 


28.5 
36.6 


46.6 
57.3 


104.1 
136.6 


ST-I-DT obs. ratio to SM 
ST-I-DT exp. ratio to SM 


10.1 
17.1 


13.4 
22.1 


22.6 
32.0 


38.6 
51.0 


93.4 
116.7 



TABLE 11: Observed and expected 95% C.L. limits 
on the cross section times branching fraction a x B, 
where B = B{H — > bb) and a is in pb, for different 
Higgs boson mass values, for single and double ^-tagged 
events, and ST-I-DT combination in the WH — > ii/bb 
channel, with £ = e or ^. The corresponding ratios to 
the predicted SM Higgs production cross section are 
also given. 
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account the difference in integrated luminosity. 

With the hmits from the WH channels reported 
above, we now turn to the combination of these 
with limits previously obtained from other chan- 
nels. We combine our new WH results with all 
the other direct searches for SM Higgs bosons pub- 
lished by D0. These are searches for Higgs bosons 
produced in association with vector bosons {pp — *■ 
ZH -> lyDbb/Ubb 4, 5], pp ^ WH WWW 0) 
or singly through gluon-gluon fusion {pp — s- _ff — s- 
WW |7[). The searches were conducted with data 
collected during the period 2003-2005 and cor- 
respond to integrated luminosities ranging from 
0.30 fb-i to 0.45 fb-i. They are separated into 
twelve final states (adding to the four WH final 
states combined earlier) and referred to as analyses 
in the following. Each analysis is designed to iso- 
late a particular final state defined by a Higgs bo- 
son production and decay mode. To ensure proper 
combination of signals, the analyses were designed 
to be mutually exclusive. 

The sixteen analyses are categorized by their 
production processes and outlined in Table IIHI 
When possible, we search for both H ^ bb and 
H WW decays. For the H ~* bb decays, 
we conduct separate ST and DT analyses, ex- 
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FIG. 4: 95% C.L. cross section upper limit (and cor- 
responding expected limit) on (j{pp WH) x B{H 
bb) (W boson decaying into a lepton -I- neutrino and 
Higgs boson into bb) vs. Higgs boson mass, compared 
to the SM expectation. The published DO e channel 
observed results, based on an integrated luminosity of 
0.17 fb~^ and the CDF {e + fj, channels) results with 
0.32 fb~^ and 1.0 fb"^ are also shown. 



Channel 


L 


Final 


Ref. 






(fb"^) 


Variable 




WH 


eubb, ST/DT 


0.43 


Dijet mass 




WH 


^ fiiybb, ST/DT 


0.45 


Dijet mass 




WH 


iubb, ST/DT 


0.30 


Dijet mass 


[41 


ZH -* vubb, ST/DT 


0.30 


Dijet mass 


[41 


ZH ' 


^ifibb, DT 


0.37 


Dijet mass 


m 


ZH eebb, DT 


0.45 


Dijet mass 


m 


WH 


^ lVWlV(e±e±) 


0.45 


LH discriminant 


m 


WH 


WWW{e^H^) 


0.43 


LH discriminant 


[6] 


WH 


WWWifi^fi^) 


0.42 


LH discriminant 


m 


H 


WW (ee) 


0.33 


Aip{e,e) 


\7] 


H ^ 


WW (e^i) 


0.32 


A<p{e,fi) 


\7] 


H 


WW (/i^i) 


0.30 


A(f{n,fi) 


17] 



TABLE III: List of analysis channels, corresponding in- 
tegrated luminosities (L), final variables for the search, 
and references. LH stands for likelihood. 



cept for ZH — ^ £~^i~bb analyses where only the 
DT analysis has been performed. The decays of 
the vector bosons further define the analyzed final 
states: WH evbb, WH pLvbb, ZH eebb, 
ZH iifibb, and ZH — > vDbb. There is a sizeable 
amount of WH tvbb signal that can mimic the 
ZH — > uDbb final state when the lepton is unde- 
tected, or when the lepton is a r decaying hadron- 
ically. This case is treated as a separate WH anal- 
ysis, referred to as WH ^vbb. 

We also include the analysis of WH 
WWW final states when the associated W bo- 
son and the same-charged W boson from the 
Higgs boson decay leptonically, thus defining six 
final states: WH We'^ve^v, We^vii^v, and 



Wfi^iyfi^i', which are then grouped into three 
analyses: e"^ e^ , fj.^ , and e'^fi'^. All decays of 
the third W boson are included. 

In the case of ^ H ^ WW production, 
we again search for leptonic W boson decays with 
three final states, WW evev, ev^v, and ^v^v. 
For the gluon-gluon fusion process, H ^ bb de- 
cays are not considered due to the large multijet 
background. 

As before, we combine results using the CL^ 
method. Systematic uncertainties are treated as 
uncertainties on the expected numbers of signal 
and background events, not on the outcomes of 
the limit calculations. This approach ensures 
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that the uncertainties and their correlations are 
propagated to the outcome with their proper 
weights. The method used here utiHzes binned 
final-variable distributions rather than a single- 
bin (fully-integrated) value. In the case of the 
H ^ bb analyses, the final variable used for limit 
setting is the invariant dijet mass, as shown for 
the WH channel in Fig. [2l In the case where 
H WW, the Higgs mass cannot be directly 
reconstructed due to the neutrinos in the final 
state. Thus, the WH WWW analysis uses 
a likelihood (LH) discriminant formed from topo- 
logical variables as a final variable while the 
pp ^ H ^ WW analysis uses the separation in (p 
between the final state leptons Ai^(£i, £2) 0- Each 
signal and background final variable is smoothed 
via Gaussian kernel estimation jl^l . 

Both signal and background systematic uncer- 
tainties vary for the different analyses. Here we 
summarize only the largest contributions, referring 
to the original publications for details. All analyses 
carry an uncertainty on the integrated luminosity 
of 6.1%. The H ^ bb analyses have an uncertainty 
on the 6-tagging rate of (5-7)% per tagged jet. 
These analyses also have an uncertainty on the jet 
energy calibration and acceptances of 8-10%. For 
the H WW and WH WWW analyses, the 
largest experimental uncertainties are associated 
with lepton measurement and acceptances. These 
values range from (3-8)% depending on the final 
state. The largest contribution for all analyses is 
the uncertainty on the background cross sections 
at (6-19)% depending on the background. The 
uncertainty on the expected multijet background 
is dominated by the statistics of the data sample 
from which it is estimated, hence is uncorrelated 
between analyses. The systematic uncertainties for 
the background rates are generally several times 
larger than the signal expectation itself and are 
thus an important factor in the calculation of lim- 
its. As such, each systematic uncertainty is folded 
into the signal and background expectations via 
Gaussian distribution. Correlations between sys- 
tematic sources are carried through in the calcula- 
tion. All systematic uncertainties originating from 
a common source, see Table IIV[ are taken to be 
correlated. 

To minimize the effect of systematic uncertain- 
ties on the search sensitivity, the individual back- 
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6-tagging (%) 


10(7) 


12 
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Background cr (%) 


6-19 
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TABLE IV: List of leading correlated systematic un- 
certainties. The values for the systematic uncertain- 
ties are the same for the ZH — > i/i/bb and WH 
^vbb channels. Each uncertainty is considered to be 
100% correlated across channels. The correlated sys- 
tematic uncertainty on the background cross section 
(cr) is itself subdivided according to the different back- 
ground processes in each analysis. 

ground contributions are fitted to the data ob- 
servation by minimizing a profile likelihood func- 
tion (ill . The fit computes the optimal central val- 
ues for the systematic uncertainties, while account- 
ing for departures from the nominal predictions by 
including a term in the function which sums the 
squared deviation of each systematic uncertainty 
in units normalized by its ±lcr uncertainties. A 
fit is performed to the background-only hypothesis 
and is constrained to bins with a signal expectation 
smaller than 4% of the total expected background. 

To set limits on Higgs boson production (cr x 
B{H X)) the sixteen analyses are first grouped 
by final state to produce individual results. We 
then group channels by production modes to form 
combined results and study their respective sen- 
sitivities. The individual analyses are grouped to 
form the LLR distributions shown in Fig. [5] for 
(a) all WH searches, with H ^ bb (ST, DT) in 
the low mass range {uih = 105 — 145 GeV), (b) 
all ZH searches (ST, DT) in the same low mass 
range, (c) all WH WWW searches, over an 
extended mass range {niH = 120 — 200 GeV), and 
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FIG. 5: LLR distributions obtained with the CLs 
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for the direct production channel, H — > WW . See text 
for details. 
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FIG. 6; LLR distributions obtained with the CLs 
method for the combination of all channels. See text 
for details. 



(d) all H WW searches, over the full mass range 
{niH — 100 — 200 GeV). We then combine groups 
(a)-(d) over the full mass range, as shown in Fig.[6l 
We also compute our results in terms of the ratio 
of the limits to the SM cross section axB{H X) 
as a function of Higgs boson mass. The SM pre- 
diction for Higgs boson production would therefore 
be excluded at 95% C.L. when this limit ratio falls 
below unity. Table |V] shows the expected and 
observed 95% C.L. cross section limits and their 
ratios to the SM for the WH and ZH analyses in 
the mass range niH = 105 — 145 GeV. Table IVII 
shows the same information for WH WWW 



Higgs mass [GeV] 


105 


115 


125 


135 


145 


WH observed a x B 
WH expected a x B 


1.60 
2.83 


1.49 
2.38 


1.57 
2.22 


1.56 
1.89 


1.65 
2.17 


ZH observed a x B 
ZH expected a x B 


2.41 
2.21 


2.23 
2.02 


1.97 
1.73 


1.77 
1.52 


3.21 
2.65 


WH observed ratio to SM 8.4 11.7 19.8 36.7 87.2 
WH expected ratio to SM 14.9 18.6 28.1 44.5 114.7 


ZH observed ratio to SM 
ZH expected ratio to SM 


21.1 
19.4 


28.5 
25.9 


40.0 
35.2 


66.0 
56.6 


263.6 
217.4 



TABLE V: Observed and expected 95% C.L. upper 
limits on the cross section times branching fraction 
a X B, where B — B{H — » bb), and a is in pb, for dif- 
ferent Higgs boson mass values, for the WH and ZH 
combined channels (WH includes the leptonic chan- 
nels, and the case where the charged lepton is not de- 
tected; ZH includes the ee,/i/x, and ui/ channels). 
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Higgs mass [GeV] 


100 


110 


115 


120 


130 


140 160 


180 


200 


WH WWW observed ax B 








11.27 4.41 


1.57 0.09 


0.010 


0.004 


WH WWW expected ax B 


- 


- 


- 


10.78 


3.53 


1.30 0.07 0.007 0.003 


H — > WW observed a x B 


10.79 


5.61 




6.07 


5.94 4.24 3.69 


4.07 


3.25 


H WW expected a x B 


8.94 


6.31 




7.74 


6.18 


5.25 3.58 


3.40 


3.98 


WH WWW observed ratio to SM 








110.7 


74.7 


53.7 46.1 


62.1 


89.6 


WH WWW expected ratio to SM 








105.9 


59.8 


44.7 34.4 


44.6 


60.8 


H WW observed ratio to SM 


636.4 


98.9 




46.1 


26.4 


14.0 9.9 


15.4 


22.2 


H WW expected ratio to SM 


527.5 


111.2 




58.8 


27.4 


17.3 9.6 


12.8 


27.2 


DO observed ratio to SM 


5.5 


7.1 


8.5 


10.5 


14.2 


12.8 10.2 


16.1 


23.7 


DO expected ratio to SM 


8.7 


10.8 


12.1 


14.3 


15.7 


13.8 9.0 


12.1 


23.5 



TABLE VI: Observed and expected 95% C.L. upper limits on the cross section times branching fraction a x B, 
where B — B{H WW) and a is in pb, for different Higgs boson mass values, for WH — > WWW and 
H WW . The ratios to the predicted values of the SM Higgs production cross section for these channels and 
for the full DO combination, are also given. 



and H WW over the full mass range. The ra- 
tios to the SM obtained with the full combination 
are also given and show the gain obtained by using 
the full information, compared to the individual 
channels. 

The expected limits for the cross section times 
branching fraction for the four groups of analyses 
(a)-(d) and for the full combination, relative to the 
SM expectations, are shown in Fig.[71 For the full 
combination of all analyses, the expected and ob- 
served cross section times branching ratio, relative 



to those for the SM, are shown in Fig. [5] Com- 
pared to an earlier simulation study of the Higgs 
boson search sensitivity conducted prior to Teva- 
tron Run II [24| . our current analyses have added 
new channels, have extended the mass range, and 
show a more uniform sensitivity for 110 < mn < 
190 GeV. 

In summary, we have presented new 95% C.L. 
limits on the WH e / fivhb production cross sec- 
tion times branching fraction which range from 1.6 
to 1.9 pb for 105 < tuh < 145 GeV. For compar- 




FIG. 7; Ratios of the expected limit on the Higgs bo- 
son production cross section times branching fraction 
to the SM expectation, for the different channel groups 
and for the full DO combination. 



FIG. 8: Ratios of the expected and observed limit on 
the Higgs boson production cross section times branch- 
ing fraction to the SM expectation, for the full DO com- 
bination. 
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ison, the expected SM cross section for uih =115 
GeV is 0.13 pb. 

We have then combined these results with aU 
previously published Higgs boson searches by the 
DO collaboration obtained with a similar luminos- 
ity (between 0.30 and 0.45 fb~^) to form new lim- 
its more sensitive than each individual limit. The 
combined observed (expected) 95% C.L. limit ra- 
tios to SM cross sections for pp WH, H bb 
range from 11.7 (18.6) at mn = 115 GeV to 
36.7 (44.5) at tuh = 135 GeV. The combined 
observed (expected) 95% C.L. limit ratios to SM 
cross sections for pp — + ZH, H ~* hh range from 
28.5 (25.9) at mn = 115 GeV to 66.0 (56.7) at 
■niH — 135 GeV. The fully combined observed (ex- 
pected) 95% C.L. limit ratio to the SM cross sec- 
tions are 8.5 (12.1) at tuh = 115 GeV, 10.2 (9.0) at 
TUH = 160 GeV, and 20.7 (16.0) at tuh = 190 GeV. 

These limits and ratios will decrease in the near 
future with the additional luminosity recorded at 
the Tevatron; more than 2 fb~^ are currently be- 
ing analyzed. New techniques are being developed 
to improve the sensitivity through advanced multi- 
variate techniques, neural-network 6-tagging, and 
improved di-jet mass resolution. In addition, an 
anticipated combination with the results from the 
CDF collaboration would yield an increase in sen- 
sitivity of about 40%. With the total expected in- 
tegrated luminosity (6-8 fb"""^), the Tevatron is ex- 
pected to provide sensitivity to the standard model 
Higgs boson beyond the current LEP limit [2^ 
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